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Abstract 

We consider CP asymmetries in neutral B meson decays to rj'Kg, (j>Ks, and K + K~ K$- We 

Oh! 

use SU(3) relations to estimate or bound the contributions to these amplitudes proportional to 

D '. 

V^Kis- Such contributions induce a deviation of the Sf terms measured in these time dependent 
' CP asymmetries from that measured for ipKs- For the K + K~ K$ mode, we estimate the deviation 

to be of order 0.1. For the n K$ mode, we obtain an upper bound on this deviation of order 0.35. 
For the 4>Ks mode, we have to add a mild dynamical assumption to the SU(3) analysis due to 
insufficient available data, yielding an upper bound of order 0.25. These bounds may improve 
significantly with future data. While they are large at present compared to the usually assumed 
Standard Model contribution, they are obtained with minimal assumptions and hence provide more 
rigorous tests for new physics. If measurements yield \Sf — S^k\ that are much larger than our 
bounds, it would make a convincing case for new physics. 
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I. INTRODUCTION 



Recent measurements of CP asymmetries in neutral B meson decays into final CP eigen- 
states test the Kobayashi-Maskawa mechanism and probe new sources of CP violation. The 
time dependent asymmetries depend on two parameters, Sf and Cf (/ denotes here a final 
CP eigenstate): 

• 4,(i) = r( B g„ ys (^/) + r( B g hy , (i )^/) = -C/cos(A mfli ) + S/Sm( A mB( ) . (i) 

CP violation in decay induces C/, while CP violation in the interference of decays with and 
without mixing induces Sf. (The contribution from CP violation in mixing is at or below 
the percent level and can be safely neglected with the present experimental accuracy.) 

If the decay is dominated by a single weak phase, Cf ~ and the value of Sf can be 
cleanly interpreted in terms of CP violating parameters of the Lagrangian. This is the case 
for decays which are dominated by the tree b — > ccs transition or by the gluonic penguin 
b — > sss transition. If one neglects the subdominant amplitudes with a different weak phase, 
the CP asymmetries in these two classes of decays are given by Sf = —77/ sin 2/3, where 
rjf = +1(— 1) for final CP-even (-odd) states and (3 is one of the angles of the unitarity 
triangle. In particular, in this approximation, the CP asymmetries in the two classes are 
equal to each other, for example, S^k = S^k- A strong violation of such a relation would 
indicate new physics [1]. 

Our aim in this paper is to quantify this statement with minimal assumptions for three 
modes of interest: 4>K S) r)'K s and K + K~K S - We would like to estimate or to find bounds 
on the deviations of the corresponding asymmetries from S^k s thai are (hadronic-)model 
independent. The ingredients of our analysis are SU(3) relations and experimental informa- 
tion on related modes. We will be able to carry out this program to the end for S v >k s - As 
concerns S^Ks, we derive SU(3) relations that can, in principle, lead to model independent 
bounds. In practice, however, some experimental information is still missing. Nevertheless, 
by using a mild dynamical assumption, we obtain a bound for this mode too. The situation 
is more complicated for Skkk, where we point out some subtleties in the interpretation of 
the experimental results. For this mode, however, we are able to estimate (rather than just 
bound) the deviation of the extracted asymmetry from sin 2(5 in the Standard Model by 
using U-spin relations and experimental data. 
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II. EXPERIMENTAL AND THEORETICAL BACKGROUND 



The CP asymmetry in B — > ipK s decays (and other, related, modes that proceed via 
b — > ccs) has been measured, with a world average [2] of 

S^ Ks = +0.734 ±0.054 [3, 4], 

CVs = +0.05 ±0.04 [3, 4]. (2) 

The value of S$k s is consistent with predictions made on the basis of other measurements 
of the CKM parameters (Am B , Am Bs , e K and tree level decays). 

CP asymmetries have also been searched for in three modes that are dominated by b — > 
sss gluonic penguin transitions: 

S lfKs = +0.33 ±0.34 [5, 6], 

C lfKs = -0.08 ±0.18 [5, 6], (3) 
S+ks = -0.39 ±0.41 [5, 7], 

C+ks = ±0.56 ±0.44 [5], (4) 

Sk+k-Ks = ±0.49 ±0.44^ [5], 
C K+K - Ks = ±0.40 ±0.341°;*; [5]. (5) 

The Standard Model predicts that in these modes —rjfSf = S^ Ks and Cf = to a good 
approximation. The statistical errors in Eqs. (3)-(5) are too large to make any firm conclu- 
sions. It is clear, however, that there is still much room left for deviations from the Standard 
Model generated by possible new physics in b — > s transitions. 

The Standard Model amplitude for these three decay modes can be written as follows: 

Af = A(B° - /) = V: b V cs a) ± V: b V us a] . (6) 

The second term is CKM-suppressed compared to the first one since 

'V* h V, 



V v; b v cs , 



v* b v us 



V cb V cs 



sin 7 = 0(A 2 ), (7) 



where A = 0.22 is the Wolfenstein parameter. It is convenient to define 



t - V ub v us u, f 

* f = V*V n c ' { } 
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and thus rewrite the amplitude of Eq. (6), 

A f = VZV es a}(l + t f ). (9) 

The SU(3) analysis we carry out allows us to bound To first order in this quantity, the 
deviation of the asymmetry from sin 2/3 is given by [8, 9] 

—rjfSf — sin 2(3 = 2 cos 2/5 sin 7 cos 5/ , (10) 

where 5/ = arg{af/a C f). The £/ parameter characterizes also the size of Cf. 

C> = -28^7^ 101, (11) 

Note that Sf can be determined from tan 5/ = {r/fSf+sm 2(3) /{Cf cos 2(3), while the following 
(^/-independent) relation between Sf, Cf and may become useful in the future: 

Cf + [{rifS f + sm2[3)/co S 2[3} 2 = 4 sin 2 7 l^] 2 . (12) 

The crucial question, when thinking of the deviation of —rjfSf from sin 2(3, is the size of 
a U f/a C f. While a C f is dominated by the contribution of b — > sss gluonic penguin diagrams, 
a U f gets contributions from both penguin diagrams and b — > uus tree diagrams. For the 
penguin contributions, it is clear that \aj/a C f \ ~ 1. (The a C f term comes from the charm 
penguin minus the top penguin, while the up penguin minus the top penguin contributes to 
a j.) Thus our main concern is the possibility that the tree contributions might yield \a U f/a C f\ 
significantly larger than one. 

For final states with zero strangeness, /', we write the amplitudes as 

Af , = A{B ° -> f) = v; b V cd b c f , + V: b V ud b} . (13) 

Here neither term is CKM suppressed compared to the other. We use SU(3) flavor symmetry 
to relate the a u f ' c amplitudes to sums of b U ff . While similar SU(3) relationships have been 
explored elsewhere [10-14], most of our results and applications are new. 

The SU(3) relations, together with the measurements or upper bounds on the rates 
for the non-strange channels plus the measured rate for the channel of interest yield an 
upper bound on \^f\. Let us first provide an intuitive explanation of this. The decays to 
final strange states, /, are dominated by the a C f terms. Those to final states with zero 
strangeness, /', are dominated by the bj, terms. Thus we can estimate \a f \ and \b U f,\ from 
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the measured branching ratios (or the upper bounds on them). Then the SU(3) relations 
give upper bounds on certain sums of the bf, and amplitudes from the extracted values 
of a j and b^,, respectively. This then gives a bound on \df/a C f\, and consequently on \£f\. 
We can also check the self-consistency of the analysis, namely that \aj\ < \Af\/\V ub V us \ and 
\b C f,\ < \ Afi\/\V c bV c d\- However, as we show below, the assumptions made in this paragraph 
can be avoided entirely. 

The SU(3) relations actually provide an upper bound on \V* b V cd a c f + V* b V ud aJ\, in terms 
of the measured branching ratios of some zero strangeness final states (or limits on them). 
Therefore, without any approximations, we can bound 



if 



V us V: b V cd a) + V* b V ud a] 
x 1 1 



if + (V us V cd )/(V ud V cl 



(14) 



V ud V* b V cs a) + V: b V us a} 

If the bound on is less than unity, then it gives a bound on We work to first 

order in if, since the naive expectation is = 0(\ 2 ). At the present state of the data the 
bounds we obtain on if are significantly larger that A 2 , so we also work in the approximation 
A 2 <C £f < 1. This is appropriate because we want to constrain the possibility \a U f/a°f\ ^> 1. 
Therefore, we take ~ £/ in what follows (although this approximation should not be 
made once the bounds on are of order A 2 ). 

The SU(3) decomposition of aj and bj, is identical with that of and b c fl although the 
values of the reduced matrix elements are independent for the u- and the c-terms. The 
SU(3) decomposition is given in Appendix A for the channels discussed in this paper. We 
use the notation a(f) = a U j C and b(f') = b U jf for equations that apply for both cases, with 
either all u or all c upper indices. Our normalization of the various amplitudes is the same 
as that of Ref. [11]. It corresponds to Y — \A\ 2 independent of whether the final particles 
are identical or not. 

The contributions to and b c j, come from penguin diagrams or the tree b — * ccq transition 
plus some form of rescattering (such as D-exchange) to replace the cc with lighter quark 
flavors. Aside from small electroweak penguin contributions, there is only an SU(3) triplet 
term in the Hamiltonian for these amplitudes. Neglecting electroweak penguins would result 
in additional SU(3) relations between the a c f and b C f, terms. We do not make such an 
approximation in our analysis, but it might be useful for other purposes. 
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III. THE CP ASYMMETRY IN B — > rj'Ks 



A. SU(3) relations 



The CP asymmetry in B — > rj'K s is expected to yield a less accurate measurement of 
sin 2(3 than the ipK s mode. The reason is that, while this decay is dominated by gluonic 
penguins, there are CKM-suppressed tree contributions that induce a deviation from the 
leading result. Nevertheless, it was argued in Ref. [15] that this deviation is below the two 
percent level. The argument was based on relating the tree contributions in B — > rj'K 
and B — > nn decays. While this may be a reasonable hypothesis, it is based on neither 
approximate symmetry nor obvious dynamical assumptions. In this section we derive a 
more rigorous (though weaker) bound on the "problematic" subleading contribution. 

The results of Appendix A imply the following amplitude relations: 



a(viK°) = -^=b( Vl 7r°) + ] jh( VlV s), 
a( Vs K°) = -^=% 8 7r°) - ^ [&(7rV) - b( mm )} . 



(15) 



2y/2 v '°" ' 4 

The states 771 and i] 8 transform as a singlet and an octet of SU(3), respectively. They are 
related to the physical 7/ and 77 states through an orthogonal rotation: 



77 = c?7i - si]s, 



77 = S?7i + C 77 8 , 



(16) 



where s = sin9 m > and c = cos9 m >. Most extractions of the mixing angle, 6 vrj >, vary in the 
10 — 20° range, and we will use m r = 20° in our numerical calculations [16]. 
In terms of physical states, we obtain from Eq. (15) the relation 



a {rfK°) = " 2 ^ (t/tt ) 



2V2 
V3s(s 2 + Ac 2 ) 



■b(vW) + 



3V3sc 2 



6(7777) 



^3c(2c 2 - s 2 



-Krrf). (17) 



4 ' ' 4 " r " 2V2 

The SU(3) analysis gives many more relations, involving both charged and neutral B decay 
amplitudes. The most general such relation, involving up to thirteen amplitudes on the right 
hand side, is given in Appendix B. With current data, Eq. (17) gives the strongest bound: 



'K s \ 



< 



v u 



V, 



ud 



0.59, 



+ 0.33, 



+ 0.14, 



\ B(r}'K°) \ B(rfK°) \ B(r}'K°) 



/3(7T°7r ) 



+ 0.53, 



B(rfrf) 



+ 0.38, 



B( m ) 



+ 0.96, 



\| B(r]'K°) ' \ B(r}'K°) ' \| B(rfK°) 



B(r)r)') 



(18) 



This bound is obtained from Eq. (17) by taking all amplitudes to interfere constructively, 
and using Q m > = 20°. The experimental upper bounds on the relevant branching ratios are 
collected in Appendix C. Using these values, we obtain 



\^ Kg \ <0.36. 



(19) 



So far only upper limits are available for many of the rates that enter in Eq. (18). Hence 
this bound is probably a significant overestimate and will improve with further data. At 
the present state of the data, we do not consider it necessary to be concerned about SU(3) 
breaking corrections. Eventually, there may be sufficient data to fix all the amplitudes a U j c , 
including their relative phases. At that point a much stronger bound can be expected, and 
allowance for SU(3) breaking corrections will need to be made. 

Using the known CKM dependence, the bound of Eq. (19) can be translated into a bound 
on the hadronic parameters, 

< 18. (20) 



l T)'K° 

This bound is much weaker than most theoretical estimates. Since the amplitudes involved 
in Eq. (17) carry different strong phases, we do not expect that they all add up coherently, 
as assumed in (18). A more plausible (though less rigorous) estimate would be that the 
left hand side of Eq. (18) is unlikely to be larger than the largest term on the right hand 
side. This estimate would give |^'^ s | < 0.14 instead of 0.36. Clearly, more data could 
significantly improve these bounds. 

The same set of SU(3) relationships can be used to carry out a similar analysis for a 
number of other modes. The relevant general relationship for i]Ks is given in Appendix B. 
Once experimental data on the asymmetry in this mode is available, it will be interesting to 
use this relationship to obtain a similar constraint on ^Ks- 

B. Using charged modes and a dynamical assumption 

One can obtain a similar bound on the ratio oif/off for the charged mode, / = i]'K + . 
The experimental situation is such that this bound is significantly stronger than the one in 
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the neutral mode. The SU(3) relations for the decompositions of aj and bj, in B + — > PP 
decays are also given in Appendix A. They lead to the following relations: 

a( Vl K + ) = %itt + ), 

a( m K+) = -L b(n + n°) - -L h (R»K + ), (21) 

and 

V6b( m n + ) = V2b(n + n ) + 2b(K°K + ). (22) 

This last equation allows us to bound $ >v 'k+ by many different combinations of decay 
modes. The most general relationship that involves only B + decay modes is 

I its+\ (3-a;)cs + (x-l)s 2 + 2c 2 , + 
a{r]K + ) = %7T + ) + %7T + ) 

+£^K* + i?) + ^b(K°K+), (23) 

where, as before, c and s parameterize 77 — 77' mixing. The parameter x is free: it allows 
us to choose, based on the state of the data, the optimal (that is, the most constraining) 
combination of amplitudes. With the branching ratios collected in Appendix C, we find that 
at present x = 3 in Eq. (23) gives the strongest constraint, 

l^jr+l < 0.09. (24) 

While the a c f amplitudes for the charged and neutral rj'K modes are the same, the a u f are 
not. The non-triplet contributions coming from the tree b — > uud terms in the Hamiltonian 
cause the differences. (This can be seen in the SU(3) relations in Table II of Appendix A.) If 
we examine the quark diagrams for these two channels we find that cfy K+ has a color-allowed 
tree diagram contribution, while a^, Ks only arises from a color-suppressed tree diagram or 
penguins. Our dynamical assumption is that the color-suppressed a%j K is not bigger than 
the color-allowed cfy K+ . This could only be violated by an accidental cancellation between 
two terms that are formally different orders in 1/N C . By making this mild assumption, we 
improve the bound on (, V 'K S by more than a factor of three over that given by the pure SU(3) 
analysis. 

There are experimental tests that could indicate that 0^,^-+ is small compared to o^^- s - 
First, if direct CP violation in the neutral mode were established to be large, \C v > Ks \ <jL 1, 
it would place a lower bound on \a%>K s / a n>K s \- Second, if the difference of the neutral and 
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charged B — > rj'K rates were sizable, given our strong bound on ^wk+I^kA-i ^ would 
imply a large \o^K s / a rfK s \ with a relative strong phase that is not close to n/2. If either 



of these measurements violate the upper bound on \a>rfK+/ a r)'K+\i ^ would suggest either an 
accidental cancellation in the charged mode, or, more interestingly, possible new physics. 
This would make it important to improve the direct neutral mode test for new physics. 

Eq. (3) shows that C v 'k s is consistent with zero and the data in Appendix C show that 
the ratio of charged and neutral rates is not necessarily much different from one. However, 
this does not validate our assumption. For example, a small strong phase, 5 v >k ~ 0, and a 
large weak phase, 7 ~ tt/2, would make direct CP violation small and induce approximately 
equal rates, independently of the size of o^, K+ / d^, Ks . 

IV. THE CP ASYMMETRY IN B — > cf>K s 
A. SU(3) relations 

A similar analysis can also be applied to B — > 4>K S . Again, the existence of CKM- 
suppressed contributions induce a deviation from the leading result, and our goal is to 
constrain that effect using SU(3) related modes. Here it is usually assumed that these 
corrections are not large, since the b — > uus tree diagram can only contribute via rescattering 
to the <j) final state which is pure ss. Thus it was generally argued that the deviation of 
S^Ks from sin 2/3 is likely to be no larger than 0(X 2 ). Ref. [12] proposed an SU(3)-based 
relation that can potentially bound this deviation, however, it involves an implicit dynamical 
assumption. In this subsection we present exact SU(3) relations that can, in principle, give 
a model independent bound. In the next subsection, we explain the dynamical assumption 
that leads to the bound of Ref. [12] and update it with current data. 

The SU(3) decomposition of aj and bj, for final states composed of a vector and a pseu- 
doscalar meson is given in Appendix A. These results imply the following relations: 




a(4> 8 K°) 




(25) 
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The states 0i and 08 transform as a singlet and an octet of SU(3), respectively. They are 
related to the physical and ui through an orthogonal rotation: 

= - ^| 08, U = ^01 + ^08, (26) 

which defines the as a pure ss state. Thus, in terms of physical states, we obtain 
a(<pK°) = l - [b(K^K°) - b(K*°K*)} + \\j\\cb{H ~ sb{H)} 

+ ^ [<*M - sb(urf)] - ^ [cb(p°rj) - sb(p°rj')} 
+\b{p^)-\b{^)-^=b{^). (27) 

This relation could give a bound on ^ Ks in a similar fashion as Eq. (18) in Section III A. 
However, a survey of the experimental data shows that currently no useful bound can be 
obtained from Eq. (27). While it is possible, using SU(3) relations, to replace some modes 
that occur on the right-hand side of Eq. (27) with a combination of others, there is no 
relation that yields a bound on £<^ s below unity at present. 

We conclude that, while in the future it will be possible to use relations such as (27) to 
constrain a^ R0 / a^ K0 in a model independent way, it is impossible to do so with current data. 



B. Using charged modes and a dynamical assumption 

For the charged mode / = (pK + , one can similarly obtain a bound on the ratio C</)K+ 
based purely on SU(3). The SU(3) relations for B + — > VP decays are given in Appendix A. 
They lead to the following relations: 

a{4>iK + ) = &(0i7r + ), 

a(0 8 K+) = 6(0 8 vr + ) - sj^b(K^K + ). (28) 
In terms of physical states we thus obtain 



a((j)K + ) = 6(07T + ) + b(K*°K + ). (29) 
Using the experimental upper limits on rates collected in Appendix C, we obtain 

|£^+|< 0.25. (30) 
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Once again there is no immediate relationship between the a^ K+ and a^ Ks . They differ 
by non-triplet Hamiltonian contributions arising from the tree-type b — > uus transition (and 
a small but similar effect from electroweak penguins). To use the above result as a bound 
on (,<f,K s requires an additional assumption that o^ Ks is not much larger than a^ K+ . Here 
we cannot readily justify this assumption, although we know of no reason why it should not 
hold. Because the is a pure ss state, there is no order N% tree contribution to a^ K +, as 
there was in the case of a^ K+ . The tree b — > uus contribution must undergo a rescattering in 
order to contribute. This brings it to be an order N c term, at the same level as the penguin 
contributions. Traditional analyses of this channel assume that the rescattering contribution 
is negligible compared to the b — > sss penguin terms, in which case the charged and neutral 
a 4>K would be approximately equal. Indeed, only an accidental cancellation could make 
|a^+| much smaller than |a$# s |. 

With this assumption, the bound for the charged mode also applies for the neutral mode. 
In this case, there is presently no pure SU(3)-based bound, and so this assumption is nec- 
essary to obtain any result. The bound in Eq. (30) was applied to the neutral mode in 
Ref. [12]. We have shown here that implicit in this bound is the assumption that there is 
no accidental cancellation of the two contributions to the charged mode. 

As in the rj'K modes, here too there is no evidence either for large C^k s or f° r large 
difference between B((j)K ) and B(4>K + ). If such evidence existed, it would indicate that 
|o£k-+| is small compared to \a^ Ks \. This could indicate that there is indeed a cancellation 
between two independent contributions to a^ K+ , or be a harbinger of new physics effects. 
Data on the non-strange neutral modes would then be desirable, as it could distinguish these 
two possibilities. 

V. THE CP ASYMMETRY IN B -> K+K~K S 

The B — > K + K~K S decay is different from the other decay modes discussed in one, very 
important, aspect: since the final state is three-body, it does not have a definite CP. Thus, 
the CP asymmetry is diluted with respect to sin 2(3. To extract the value of sin 2(3 from this 
measurement, one has to know the relative fractions of CP-even and CP-odd final states. 
The BELLE collaboration employed a beautiful isospin analysis for this purpose [17]. The 
accuracy of the isospin analysis affects the accuracy with which the true Skkk (that is, 
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Skkk for a final KKK state with a definite CP) is determined. Consequently, the relation 
between the experimental value, S^kkj and sin 2(3 is more complicated. We will first analyze 
the accuracy of determining Skkk from S^kk and then the deviation of Skkk from sin 2f3. x 

A. Isospin analysis 

The B — > K I K J K L decays (where I, J, L = {+, — , 0, S} specify the kaon states) involve 
an initial 1=1/2 state and final 1=1/2 and 3/2 states. There are five independent isospin 
amplitudes, A\ , Af , ^4|' and where the lower index denotes the isospin representation 
of the Hamiltonian and the upper index gives that of the final state. (We follow the notation 
of the previous sections, instead of using isospin labels.) For the isospin-doublet final states, 
the representation 2 denotes where the two S — — 1 mesons are in an isospin-singlet state, 
while 2' denotes where they are in an isospin-triplet. Defining 

A IJL ( Pl ,p 2 ,p 3 ) = A(B - K I (p 1 )K J (p 2 )K L (p 3 )), (31) 

we obtain the isospin decompositions given in Appendix D. 

Let us start by neglecting the tree contributions, as was done in the BELLE analysis. 
This corresponds to A\ — A% — A\ — 0. Then, the following amplitude relations arise: 

A)o+ = A + ^ , A +00 = A -+ , A 000 = A + - + . (32) 

When integrating over phase space, the contribution from the interference between A\ and 
A\ vanishes. (Thus, although Ajjl is not invariant under I <-> L, the rates Tjjl and the 
branching ratios Bjjl are.) Consequently, the equalities of the following rates are predicted: 

r+-o = r +0 o , r + _ + = r oo • (33) 

Branching ratios of four B — > KKK decays have been measured [17]: 

B + _ + = (3.30 ± 0.18 ± 0.32) x 10" 5 , 

B + _ = (2.93 ± 0.34 ± 0.41) x 10" 5 , 

B +ss = (1.34 ± 0.19 ± 0.15) x 10" 5 , 

B sss = (0.43^ ± 0.75) xlO- 5 . (34) 

1 In the original version of this section, there were errors in the isospin and U-spin decompositions of the 
relevant amplitudes, which were pointed out in Ref. [18]. We agree with their results. 
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Thus the approximation in Eq. (32) that lead to (33) is not yet tested (in particular, isospin 
symmetry implies no relation between £>+_ + and £> + _o)- Because Tj 00 include both CP-odd 
and even states for the pair of neutral kaons, the measured rates in Eq. (34) are not sufficient 
to test the relations in Eq. (33). For example, the first relation in Eq. (33) becomes 

r+- s = 2 r+sx + r+ss . (35) 

We now focus on the K + K~K° and K°K°K + modes. One can write down the effective 
Hamiltonian in terms of the meson fields. The two 1 = terms are of the form (B l K i )(K : > Kj) 
where i and j are isospin indices. We can decompose the Hamiltonian into components where 
the K 3 Kj pair is either in / = even or in I = odd angular momentum state: 



H eS oc (B'K,) [x (K^Kj) l=evm + VT^ 2 (K 3 K,) l=odd \ , (36) 

The equality of the amplitudes in Eq. (32) guarantees that x is equal for the two decay 
modes, and allows the extraction of the CP even/odd fractions in the B° — > K S K + K~ 
decay from measurements of the B + — > K + K S K S decay as follows [17]. 

Consider the K + K~ sub-system in the B° — > K + K~K° decay. Charge conjugation 
exchanges K + and K~ , and Parity exchanges them again (in the center of mass frame). Thus 
the K + K~ system has CP= +1. Then, the K + K~K$ system has CP= (—1)', where / is the 
relative angular momentum between K s and (K + K~). (It also equals the relative angular 
momentum between K + and K~ .) There is then a one-to-one correspondence between 
angular momentum and CP. In particular, x 2 in Eq. (36) gives the CP-even fraction in the 
B -> K + R-K s decay. 

Next consider the K°K° sub-system in the B + — > K + K°K° decay. Bose symmetry 
implies that I = even corresponds to a final K$Ks + K^K^ state, while / = odd corresponds 
to a K S K L state. Thus, x 2 = 2r +i 9g/r +0 o- Since B +00 has not been measured, one can use 
again Eq. (33) to arrive at the following relation [17]: 

x 2 = 2— — = 0.97 ±0.15 ±0.07. (37) 

r+-o 

We learn that the K + K~K S final state is dominantly CP-even. 

From a measured value of the CP asymmetry, S C ^ K , we can deduce the value of the CP 
asymmetry for the CP-even component, S KKK , according to S KKK = S c ^l CK /(2x 2 — 1). We 
learn that in the limit that 1=1 contributions to the Hamiltonian are neglected, we have 

QCXp 

C ^KKK ( oo\ 

41 +SS/1 +-0 - 1 
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When the higher isospin contributions are taken into account, the three amplitude equal- 
ities of Eq. (32) and the two rate equalities of Eq. (33) no longer hold. There remains a 
single amplitude relation, 

A 000 + A + _ + + A +00 + A 00+ + A + _ + A)-+ = . (39) 

The angular momentum analysis is modified by the three A 3 terms. In particular, both the 
relation T + _s(/ = even)/r+_s = r +0 o(^ = even)/T +0 o and the relation r +0 o = T+-o are 
corrected by terms of 0[{A\ + \/2A§ / (y/ZA\)\ and of 0[A$ /(y/3A% )]. At present there is 
no experimental information on the size of these corrections. 

One might worry about isospin violation in the — > KK decays, since B((f) — > K + K~) m 
49% and B(4> — > K S K L ) xs 34% should be equal in the isospin limit. (This large violation 
can be understood as arising chiefly from the phase space difference for the two channels.) 
Since B(B — > <f>K) x B(<j> — > K + K~) is between 10 - 15% of £> + _ , this could give an 
additional error of up to ~ 4% on x 2 , not a very large effect. 

Note that even if the A 3 amplitudes were negligibly small and thus the isospin analysis 
to find the CP-even fraction in the KKK state very precise, it would not imply that the 
extracted CP asymmetry is equal to sin 2/3 to the same precision. The b — > uus tree con- 
tribution also has an isosinglet component, which would not affect the isospin analysis but 
would shift Skkk from sin 2(3. In the next subsection we estimate the overall effect of that 
contribution. 

B. U-spin analysis 

In the previous subsection, we used isospin symmetry to estimate the CP-even fraction 
in the K + K~ K$ final state. Isospin symmetry relates the B° — > K + K~K° mode to the 
B + — > K + K°K° mode. In this subsection we use U-spin symmetry to estimate the overall 
effect of contributions to B — > K + K~K + that are proportional to V* b V us . U-spin relates 
certain B + — > hfhjh^ modes to each other, where hij t k = K or ir. (Since U-spin is a 
subgroup of SU(3), this analysis is just a simplified form of the analysis that we have given 
for the other channels in this paper.) 

Under U-spin, B + is a singlet, while Mj = (K + , 7r + ) is a doublet. A crucial point in our 
discussion is the U-spin transformation properties of the Hamiltonian. Both the penguin 
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amplitudes, b — > (mm + Jrf + ss)g, and the tree contributions b — > img (with q — d,s) are 
At/ =1/2. Consequently, there are two U-spin amplitudes for the charged B decays into 
three charged kaons or pions. The decomposition of the various decay amplitudes in terms 
of these two U-spin amplitudes is given in Table IV. We find the following relation: 

a(K + K~K + ) = b(n + n-7i + ). (40) 

The experimental data are [17, 19] 



1T+K-K+) = (6.6 ±3.4) x 1(T 6 , 
>K + n-7i + ) = (5.7 ±0.4) x 1(T 5 , 



BrKK = 


B{B + - 


B n KK = 


-,B(B+ 


Bk-k-k 1 


= B(B+ 


B 


= B{B^ 



1-5 



(41) 



To relate the two pairs of rates in a useful way, we make our usual approximation: we take 
the KKK rate to be dominated by the a c KKK term, and the nnn rate to be dominated by 
the a" term. Then we obtain 



TVTTTT 



KKK\ — 



V u . 



B(B+ -> 7T + 7T-7r + ) . , 

~ 0.13, (42) 



V ud \ B(B+ - K+K-K+) 

Given the size of U-spin breaking effects and the crudeness of our approximations, we esti- 
mate that the corrections to —Skkk = sin 2(5 is of the following size: 

Zk+k-Ks =0.13 ±0.06. (43) 

Additional constraint on \^kkk\ can be derived from B^kk- Our amplitude relations 
imply that, in the U-spin limit, B n KK > B n7T7T /2. Consequently, 



\£,kkk\ < 



V 



ad 



\ B(B+ -> K+K-K+) 

The above analysis does not distinguish between quasi-two-body and true three-body 
contributions to the eventual three-body rate. Indeed it includes all three-body final states, 
whether or not reached by a resonant contribution. We used here the SU(3) relationship only 
for the total rates, integrated over the entire Dalitz plots. Comparisons of more restricted 
regions of the Dalitz plots would be much more subject to SU(3) breaking corrections. 
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VI. CONCLUSIONS 



\€r,'K s \ < 



Within the Standard Model, the CP asymmetries Sf in neutral B decays to the final 
CP eigenstates 4>K S) rj'K s and (K + K~ K s )cp=-i are equal to the CKM parameter sin 2(3 
measured in B — > ipK s , to a good approximation. Furthermore, the direct CP asymmetries 
Cf in these modes are expected to be small. The goodness of this approximation is different 
between the various modes and its estimate suffers, in general, from hadronic uncertainties. 
We used SU(3) relations and experimental data (and, in some cases, a mild dynamical 
assumption) to estimate or to derive upper bounds on the deviation of the Sf from sin 2/3 
and on the size of Cf. We obtained 

0.36 SU(3), 

0.09 SU(3) + leading N c assumption, 

\£,<t>K s \ < 0-25 SU(3) + non-cancellation assumption, 

\i K +K-K s \ ~ 0.13 U-spin, (45) 

where £/ is defined in Eq. (8). The approximations and assumptions that lead to these 
results are spelled out in the corresponding sections. While our bounds for the first two 
modes are considerably weaker than estimates based on explicit calculations of the hadronic 
amplitudes, they have the advantage that they are model independent. Although SU(3) 
breaking effects could be significant, our bounds for the two-body modes are probably still 
conservative because they arise from a sum over several complex amplitudes that we assumed 
to interfere constructively. Furthermore only experimental upper bounds are available for 
many of the rates that enter these bounds. As data improves these bounds could become 
significantly stronger. Certainly, if deviations from sin 2f3 are established that are larger than 
the SU(3) bounds, the case for new physics would be convincing. Since our bounds apply 
more generally to minimal flavor violation models, the new physics would have to be beyond 
this framework. Even where our results require additional assumptions, the situation here 
is better than the usual, in that we are making assumptions about non-leading corrections 
to the af amplitudes, rather than about the full aj terms. 
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APPENDIX A: SU(3) DECOMPOSITION FOR (f^\(bu)(uq)\B) AMPLITUDES 

In this appendix we give the SU(3) decomposition of matrix elements that are relevant to 
our analysis. The operator that creates a B meson containing a b quark transforms as a 3 of 
SU(3). The AB = +1 Hamiltonian, which has the flavor structure {bqi){qjqk), transforms 
as 3x3x3 = 15 + 6 + 3 + 3. Our calculations follow closely that in Ref. [20], the only 
difference being the decomposition of the Hamiltonian that can be read off from Ref. [10] . 

TABLE I: SU(3) decomposition of a(f) and b(f') for /W = 771 P 8 - 



/(') 


08 
D 15 


08 
^6 


si 




-1 


-1 


1 




3 


1 


1 


77i7T° 

ViVs 


v/3/2 


-1/V2 

-v/3/2 


-1/V2 
1/V6 


771 7T + 


3 


1 


1 
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TABLE II: SU(3) decomposition of a(f) and b(f') for /(') = P 8 P 6 



/(') 


Ail 


^15 


^1 




^3 




'IS 1 *- 


4a/r/5 


1 A/fi 


— 1 /\/fi 


— 1 A/fi 


n 


K°TT° 


1 9a/9 


1 /a/9 
1/ v z 


— 1 A/9 
1/ V z 


— 1 /a/9 
1/ V z 


n 

u 




16/5 


-1 


1 


1 





'/o J *■ 


8\/6/5 


-a/372 


l/\/6 

/ V 


-1/a/6 

/ V ^ 


o 


K+7T 


16\/2/5 


3/a/2 


-1/V2 


1/V2 





X 7T+ 


-8/5 


3 


-1 


1 







n 
u 


0/ Vo 


1 A/ 1 ? 


-1/V3 





7T 7T 


— lov z/o 


1 / /o 
1/ V 2 


1 / /o 
1/V 2 


1/(3V2) 


V2 


(Jin Tin 

'18 '18 


0\ A 1 O 


-1 /a/9 
±/ V z 


— 1 /a/9 
±/ V z 


-l/(3>/2) V2 




14/0 


i 
1 


i 
1 


1/3 


2 


K~K + 


-2/5 


2 





-2/3 


2 




-2/5 


-3 


-1 


1/3 


2 


r] 8 ir + 


4^6/5 


x/6 




a/2/3 





TT+TT 


4^/2 
















-8/5 


3 


-i 


1 






For final states composed of an SU(3) singlet and an octet, there are three reduced matrix 
elements. The reason is that there is a unique way of making a singlet from an octet plus 
any one of the three representations of the Hamiltonian and the B operator. In Table I 
we give the decomposition of a(f) = a U j c and b(f') = b'^f for /W = r]iP 8 , where r/i is the 
SU(3)-singlet pseudoscalar, and P 8 is the SU(3)-octet pseudoscalar. The matrix elements S% 
that occur in the decomposition of and bj, are independent of those for a,f and Uj,. In our 
notation, the lower index of the reduced matrix elements denotes the SU(3) representation 
of the Hamiltonian and the upper index is that of the final state. (If electroweak penguin 
contributions were neglected, the decomposition of a c j and b c j, is given by the last column, 
corresponding to if in a triplet.) The decomposition for /W = faP$, where fa is the SU(3)- 
singlet vector-meson, is the same as that for /W = r/iP 8 , with different values of the reduced 
matrix elements, S%, and the replacement r)i — > fa. 
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Final states containing two SU(3) octets can be decomposed as 8 x 8 = 27 + 10 + 10 + 
85 + 8,4 + 1. The final state composed of PsPs is symmetric, and so it transforms as an 
element of the symmetric part, (8 x 8)5 = 27 + 8 + 1. In Table II we give the decomposition 
of a(f) and b(f) for /W = P 8 P 8 ; it contains five reduced matrix elements. When the final 
mesons are different, such as /W = P 8 V 8 , where V 8 is the SU(3)-octet vector-meson, all 
six representations appear. In Table III we give the decomposition of a(f) and b(f') for 

= P 8 V 8 , which contains ten reduced matrix elements. Again, the matrix elements 
and B@ that occur in the decomposition of aj and Vj, are independent of those for a c f and 
b°f,. (If electroweak penguins were neglected, the decomposition of a C f and bf, is given by the 
columns corresponding to if in a triplet, A\,A\, B\ s , and B\.) 

Related tables have been presented in Ref. [10], however, the (bu)(us) contributions to 
strangeness changing decays were neglected. Furthermore, in that paper, when there are 
several independent amplitudes with the Hamiltonian in a given SU(3) representation, these 
contributions are not decomposed according to the SU(3) representation of the final state. 
Ref. [11] gives the SU(3) decomposition for /W = P 8 Ps in a somewhat different notation 
from ours. They do not discuss the applications we investigate here. In Ref. [14] a nonet 
[U(3)] symmetry is assumed (which can be justified in the large N c limit). It relates matrix 
elements involving 0i and (f> 8 (and, similarly, 771 and rj 8 ) that are independent of one another 
based only on SU(3). Thus, in Ref. [14], eleven amplitudes describe B — > VP decays (with 
V in a singlet or octet and P in an octet), while we need thirteen. The use of the fewer 
number of matrix elements amounts to a dynamical assumption beyond SU(3). 



APPENDIX B: SU(3) RELATIONS FOR B° -» r)'K° AND B° -> -qK 



The most general SU(3) relation between the a(r)'K°) and 6(/')'s of charged and neutral 
B decays can be written as follows: 



a(r]'K 



+ 



s 2 -2c 2 y/3s 2 (xi-x 2 ) 

2V2 2 
a/3s ^ s(xi + x 2 + 4x 3 ) 



+ 



2V2 

3V3sc 2 3sc 2 ( Xl + X2) 
4 2~V^ 



b(rj'ir ) - 

b(7v\°) - 

b(rjrj) + 



3sc v^3sc(xi — X2) 



b( V 7l°) 



2^2 2 
\/3s(s 2 + 4c 2 ) 3s 3 (x 1 +x 2 



4 2V2 
VQc(2c 2 -s 2 ) 3cs 2 (x 1 +x 2 ) 



b(r)'r)' 
6(7777') 



sx 3 b(n + n-) - sxi b(K + K~) - sx 2 b(K°K°) - sx 4 b(K°K + ) 
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TABLE III: SU(3) decomposition of a(f) and b(f') for /(') = 7 8 P<. 
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r?27 
S 15 
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B 15 
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4/3 
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+ 



y| S cx 4 %7T + ) - ^s^b^) + (V2sx 3 - b(n°n + ). (Bl) 



Here, x±, x 2 , x% and £4 are free parameters that allow us to choose between various com- 
binations of amplitudes on the right hand side of this relation. In particular, given a set 
of experimental measurements of (or bounds on) the corresponding branching ratios, we 
can vary the Xi parameters so that we get the strongest constraint. With current data, the 
optimal choice is X\ = x 2 = x 3 = x 4 = 0, which yields Eq. (17). 

The analogous relation that will allow to bound the deviation of the CP asymmetry in 
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B — > 7]K S decay (once it is measured) from sin 2(3 is: 



v/2 



6(r/7T° 



+ scx 2 %'7r+) - ^| c 2 x 2 %7r+) + | |^c 2 ^ + 3x 3 j - s 2 ^ 6(7/77) 



3s 2 c 



c 

(v 7 ^ + V2x 3 ) btfrf) + c — (Vz - ZV2x 3 ) + s 2 V3 



Km) 



- | (>/3 + ^(4^4 + a*)) &(7r°7r°) + cx 4 b(ir + n-) + -^=(x 2 - 2x A ) b(n\+) 

+ C - (x 3 - Xl ) b{K+K~) + °-{x 3 + Xl ) b(K°K*) + cx 2 b(K*K + ) . (B2) 

We have also derived the most general SU(3) relation between the a((psK°) and the 
(sixteen) of charged and neutral B decays. The relation is quite complicated and it 
does not seem likely that it will become useful in the near future, so we do not present it 
explicitly here. 



APPENDIX C: RELEVANT BRANCHING RATIOS 

The values below are collected from Ref. [16]. 



B{r]'K ) 
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B(rj'K + ) 
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B(r)r)) 


< 


1.8 x 10" 5 , 
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B(r]ri') < 2.7 x 1(T 5 , 
B(r]'r]') < 4.7 x 10~ 5 . 



(CI) 



B(cj)K ) 
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(8.1+11) x 10 


B(<j)K + ) 
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(7.91?;°) x 10 
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1.0 x 10~ 5 , 


B(v'p°) 


< 
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3.0 x 10" 6 , 


B{K + K*°) 


< 


5.3 x 10" 6 , 


B{(f)7l + ) 


< 


1.4 x 10~ 6 . 



APPENDIX D: SU(2) DECOMPOSITION FOR (KKK\(bu)(us)\B) AMPLI- 
TUDES 

In this appendix we give the isospin decomposition of matrix elements that are relevant 
to our analysis. The operator that creates a B meson containing a b quark transforms as a 
2 of isospin-SU(2). The AB = +1 Hamiltonian, which has the flavor structure {bs)(q i q, i ) ) 
transforms as either 2x2 = 1 + 3 for qi = u, d or 1 for qi = s. For final isospin-doublet states, 
the non-primed (primed) isospin amplitudes correspond to the K 1 K L sub-system being in 
an isospin-singlet (-triplet) state. A similar form of decomposition holds for the U-spin 
amplitudes in B + — > hf (pi)h~ (p2)hf (p^) decays (where hf^ = K + ,ir + and hj = ir~,K~). 
They are given in Table IV. 
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TABLE IV: Left: Isospin decomposition of Ajjl = A(B - 
decomposition of A hihjhl = A(B hf(p 1 )hj(p 2 )h^(p 3 )). 



K I (p 1 )K J (p 2 )K L (p 3 )). Right: U-spin 
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